function f = alobj(x)
    global pena;
    N_inequ = 4;
    h_inequ = 0;
    g = constrains(x);
    for i = 1:N_inequ
        h_inequ = h_inequ+ pena * (max(0,g(i))).^2;
    end
    f = obj(x)+h_inequ;
end
 %计算当前惩罚函数值
